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ROTATIONAL -TEMPERATURE DETERMINATION IN FLOWING 


NITROGEN USING AN ELECTRON BEAM 

By D. C. Lillicrap* and Louise P. Lee 
Langley Research Center 

SUMMARY 

The measurement of rotational temperature by the electron-beam fluorescence 
technique is investigated for temperatures between 78 K and 300 K by using the 0-0 band 
of the nitrogen first negative system. Existing methods for calculating the rotational 
temperature from electron-beam spectra result in values which are higher than the true 
temperature and the error increases with density and electron-beam current. The 
errors are believed to arise through neglect of the effects of secondary electrons. 
Excitation by low-energy secondary electrons can increase the population of the sparsely 
populated rotational ground states at the expense of the abundantly populated states. This 
secondary excitation is proposed as the cause of the error in the calculated rotational 
temperature and a theoretical study determines that the density and beam current govern 
the excitation rate. Rotational -temperature measurements are made in a flowing stream 
of nitrogen at a known temperature for a range of densities and beam currents. The 
product of density and beam current correlates the present data, and least -squares cubic 
polynomials are fitted through the temperature data as a function of this correlating 
parameter. An iterative procedure based on the polynomials allows more accurate rota- 
tional temperatures and, consequently, densities to be obtained with the electron -beam 
fluorescence technique. This procedure can be applied to measurements in some high- 
speed flow regimes. 


INTRODUCTION 

As flight speeds are extended into the hypersonic regime, high recovery tempera- 
tures are encountered and can result in low -density regions within the flow field. The 
study of these regions with conventional probes may lead to unacceptably large errors. 
For instance, in high-speed low-density flows, corrections to pressure measurements 
for probe interference and orifice effects are large and difficult to apply, particularly in 
boundary layers close to the surface. Consequently, electron beams are used extensively 

*NASA-NRC Research Associate, presently affiliated with Imperial College, 

London, England. 



to study this flow regime since they have no appreciable effect on the flow, and measure- 
ments can be made to within a few mean free paths of the surface. (See ref. 1.) 
Interference -free measurements of this kind would be particularly valuable in the study 
of high Mach number turbulent boundary layers. In the low -density region close to the 
surface, it has been found (refs. 2 and 3) that the static pressure appears to increase 
significantly as the surface is approached. Wallace (ref. 4) used an electron -beam probe 
to determine mean and fluctuating densities in a turbulent boundary layer but the results 
were preliminary and temperature profiles were not measured. The potentially greater 
accuracy of interference -free electron-beam measurements of both density and tempera- 
ture could provide a better understanding of the low-density regions of turbulent boundary 
layers. 

The electron-beam fluorescence technique (ref. 5) consists of directing a collimated 
beam of moderate energy electrons through the flow field and, by use of a suitable optical 
system, of observing the fluorescence produced. This paper considers only density and 
rotational -temperature measurements in pure nitrogen. The intensity of the first nega- 
tive bands can be related to the density and, if the rotational fine structure of a band is 
resolved, a rotational temperature can be obtained from the relative line intensities. For 
most conditions of interest, the rotational temperature equals the translational tempera- 
ture and, in principle, the state of the gas can be determined at any point in the flow with- 
out appreciably disturbing the flow. (See ref. 1.) In practice, the rotational temperature 
calculated by using the equation derived by Muntz (ref. 5) is too high and the discrepancy 
increases with density and the number of spectral lines used (refs. 6 and 7). This dis- 
crepancy also leads to errors in determining the number density for conditions where 
quenching of the excited states throi^h collisions is appreciable. (See ref. 8.) 

The objective of this paper is to examine the possible causes of the high calculated 
rotational temperatures and hence remove some of the uncertainties associated with den- 
sity and rotational-temperature measurements, particularly at the upper end of the density 
range where electron-beam techniques are applicable, on the order of 10^3 molecules 
per meter^. Excitation by low-energy secondary electrons is investigated, since it could 
result in a redistribution of the molecules among the ground-state rotational levels. An 
excitation model is proposed and the parameters affecting the calculated temperature are 
identified. The equations describing this excitation model involve unknown constants, and 
an experimental calibration of the technique for temperatures between 78 K and 300 K is 
performed for a range of densities and beam currents. An iterative procedure based on 
tabulated experimental results is developed to obtain more reliable values of rotational 
temperature and consequently density, provided certain conditions are satisfied. 
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SYMBOLS 


A,B,C 

By- 

CQjCijCgjC 

c 


F 

F(K’) 

f(v) 

G 

g(v,±2)j^„ 

H 

h 


constants 

rotational constant for Vj vibrational level 

g cubic polynomial coefficients 

mean molecular velocity 

experimental constant 

speed of light 

excitation function 

defined by equation (19) 

secondary-electron velocity distribution function 
defined by equation (2) 

excitation cross section for secondary electrons 
defined by equation (7) 

Planck's constant 


I(K’,K")=I(K’)=i(k',K' 2) intensity of K' rotational line 

I^Q intensity of R branch of 0-0 band 

i beam current 

J defined by equation (8) 

j current density 


K 


rotational quantum number 



K’ 


K" 




^max 

k 


N 


ft 

Nr,! 


nJb2 


1^X2 


NjX' 




R 


% 

Q[(v,v';) 


rotational quantum number for N^B^S state 
rotational quantum number for N2X^S state 
rotational quantum number for N^X^S state 
maximum number of spectral lines used 
Boltzmann constant 
nitrogen number density 
number density of K" level 

upper electronic state of nitrogen ion Ng 

ground state of nitrogen ion 

ground state of nitrogen molecule Ng 
secondary electron number density 

relative rotational transition probability for P branch 
relative rotational transition probability for R branch 
beam diameter or field-of-view width 
rotational state sum 

Franck-Condon factor for transition between V and Vj 
of two electronic states 


R(T) = 




net rate of secondary -electron excitation of K" level 


vibrational levels 
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r rate of excitation of K" level by one transition 

T temperature 

Tj^ rotational temperature 

wall temperature (equals static temperature) 

+ 2 

V vibrational quantum number for NgB S state 

Vj vibrational quantum number for state 

V electron velocity 

Vq reference electron velocity 

X(K') defined by equation (17) 

A(K) change in rotational quantum number 

A(K') change in upper state rotational quantum number 

A(K") change in ground-state rotational quantum number 

AI(K') difference between experimental and theoretical rotational line intensities 
6 rarefaction parameter 

Q defined by equation (9) 

\g secondary-electron mean free path 

V wave number of emitted radiation 

# rate of excitation defined by equation (5) 

(f) = Byn 


h^k 
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Subscripts: 


calc calculated values defined by equation (11) 

exp experimental values 

M Maxwellian distribution 

max maximum value 

min minimum value 

1,2 iteration approximations 

THEORETICAL EXCITATION MODELS 


Previous Models 

To obtain the rotational temperature, Muntz (ref. 5) used a theoretical model to 
relate the measured rotational line intensities for the first negative bands to the 
undisturbed ground-state rotational distribution. He argued that primary electrons obeyed 
the optical dipole selection rules (AK = ±1) in exciting the ground-state molecules directly 
to the state. This state decays spontaneously to the state and gives 

rise to the first negative system. If it is assumed that no other processes significantly 
affect the relative rotational line intensities, the excitation model leads to the equation 

^ -K’(K’+1)0/Tr 

I(K',K") = Constant :^%'Ge ' (1) 


where 


-2(K’+1)0/T„ 2K'0/t„ 

g _ (K' + l)e ' ^ + K'e ' ^ 


2K’ + 1 


( 2 ) 


An iterative procedure is used to solve these equations for Tj^, but the values 
obtained are too high. Robben and Talbot (ref. 9) first derived a correction curve for 
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(Tr^ as a function of the temperature but subsequently, it was shown that the error 
' 'calc 

in (Tpj) also increased with density. (See refs. 6 and 7.) Secondary electrons, 

'' “^calc 

produced when the primary electrons ionized ground-state molecules, were suggested 
(refs. 10 and 11) as the cause of the high calculated temperatures. Some secondaries 

1 n 

have sufficient energy to excite ground-state molecules to the state, but for 

these low electron energies, transitions for which AK = ±1, ±3, ±5, . . . can occur. This 

condition results in high calculated temperatures, but order-of-magnitude calculations 

show that there are insufficient secondaries with the required energy to account for the 

large errors in (T-r) . Furthermore, any appreciable direct excitation of the 

^ ^^calc 

+ 2 

N 2 B S states by secondaries would result in a nonlinear variation of band intensity 
with density, even when quenching is negligible, as can be seen from reference 10. For 
low densities, less than 1022 molecules per meter^^ the band intensity varies linearly 
with density, and Maguire (ref. 7) used this variation to argue against secondary electron 
excitation in favor of preferential quenching of the rotational states. For this model to 
work, the quenching cross sections for the rotational states have to be a complicated 
function of the temperature and some states need to be strongly quenched even for den- 
sities less than 1Q22 molecules per meter^. Yet another approach to this problem is 
the double excitation model proposed by Hickman (ref. 12) which allows transitions for 
AK = ±1 and ±3. Although the calculated temperatures are lower, this model offers no 
real advantage over Muntz's model. 

Present Excitation Model 

When a beam of moderate energy electrons passes through nitrogen, some of the 
ground-state molecules are ionized. The ions and secondary electrons diffuse out of the 
beam region, but the electrostatic field set up aroimd the beam retards the electrons and 
results in high concentrations of secondary electrons. Beinkowski and Harbour (ref. 13), 
using a formalism similar to that of Self and Ewald (ref. 14) for a discharge column, 
calculated values several orders of magnitude larger than those derived previously with 
the electrostatic forces neglected. For a 10-3 meter beam radius and a 10”^-ampere 
beam current in a gas at a density of 3.2 x 1022 molecules per meter3, the estimated 
secondary-electron concentration was 4 X lol"^ per meter3. Harbour, Bienkowski, 
and Smith (ref. 15) subsequently measured concentrations close to the predicted values 
and an electron temperature of about 1 eV. For electrons of this energy, the cross sec- 
tion for excitation of ground-state rotational levels with a selection rule AK" = ±2 is 
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approximately 10 meter2. (See ref. 16.) In view of the large secondary-electron 
concentrations and excitation cross sections, the effect on the ground-state rotational 
distribution function must be considered. The population of a particular ground state K" 
is influenced by the net result of four transitions between it and the levels K" = ±2, as 
illustrated in figure 1. These secondary transitions will cause a redistribution of the 
molecules among the ground-state rotational levels; thus, the population of the sparsely 
populated levels is increased at the expense of the more abundantly populated levels. 

This condition produces the effects observed by Maguire (ref. 7) without the need to 
assume that the abundantly populated states are preferentially quenched. The 

redistribution of the molecules will not appreciably change the total intensity of the first 
negative bands, and at low densities (less than 1Q22 molecules per meter 3 ) ^ this model 
predicts a linear variation of intensity with density, as observed experimentally. There- 
fore, the objections to earlier modifications of Muntz's model (ref. 5) do not apply to 
this model, but further comparison with experiment requires a detailed study of the effect 
of the secondary electron excitation on the relative line intensities. 

Figure 1 illustrates the complete excitation model by showing the transitions 
affecting the seventh rotational line. The K* = 7 level is populated by primary elec- 
trons from the K" = 8 and K" = 6 levels. Each of these levels is influenced by tran- 
sitions excited by secondary electrons between them and adjacent even-numbered lines. 



Figure 1.- Excitation model showing transitions for l(K' = 7). 
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To include the secondary -electron excitation, the perturbation of a particular K" level 
is assumed to be proportional to the net rate of excitation of that level Rj^n. Then, 

where C is a constant and the subscript M denotes the population for a Maxwellian 
distribution. 

Order-of-magnitude calculations based on equation (3) predict an overpopulation of 
the higher rotational states compared with a Maxwellian distribution. Calculations are 
performed for unit volume occupied by a beam of diameter p, with the assumption that 
molecules enter and leave the beam only by thermal diffusion through the edge of the 
beam. The difference between the rate of diffusion out of the beam for a particular 
K" level, and the rate of diffusion into the beam for the same level is the rate of excita- 
tion of that level per unit volume per second by the secondary electrons. It can be shown 
that 




from which 





M 


+ 1 




( 4 ) 


The ejqjression for R^, is developed later and depends on the density, secondary elec- 
tron concentration, secondary electron energy, and corresponding excitation cross sec- 
tions. Cross sections for 1 eV electrons were estimated to be two orders of magnitude 
greater than the values listed in reference 16 for 0.14 eV and are of the order of 
10"^^ meter^. By using a secondary electron concentration ng of 4 X 10^'^ per meter^ 
from reference 13, equation (4) predicts a considerable overpopulation of the high rota- 
tional states. For example, at 78 K an overpopulation of approximately 30 percent is 
predicted for the eleventh rotational level. These order-of-magnitude calculations 
demonstrate that a more complete study of the proposed secondary-excitation mechanism 
is required. 

To study the effect of secondary electrons in detail, the excitation by low-energy 
secondaries must be included in the equation for the relative rotational line intensities. 
Based on the Muntz method (ref. 5), the rate of excitation by primary electrons from 
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■+-d2'^ 


a V’^ vibrational band to a K' level of the V band of NgB S is 


$(V’,K’) = FNi ^ 


1 ^K’+l^P 

q(v,v;)| 

] Z (^K’+l^P ^K*-1 ^r) 


J 


( 5 ) 


per second per unit volume. In equation (5), F is the excitation function for an electron 
defined so that if all molecules are in the ground state, then FNi molecules are excited 
to the states per second per imit volume. (In the previously used notation 

(ref. 8), F was defined so that Fi molecules were excited to the NgB^S states per 
second per xmit volume. As a result, the so-called constants in the equations of ref. 8 
involve the density.) For a Maxwellian distribution 

2K" ^ 1 -K"(K"-H)^./Tr 


''i Qr(Tjj) 


V, 


where is the rotational state sum for the V'^ level. Therefore, from 


equation (3) 


”k" ■■ ”v,' 


2K.1 -K"(K’V1W/T„ 


Qr('^r)^„ 


CR^„ 


Substituting for and j in equation (5) and rearranging the equation yields 


$(V’,K’)=^P^jq(V’,V’i) 


v; 


N, 


«r(Tr) 


H + J 


( 6 ) 




per second per unit volume. In equation (6) 

.[(K’+l)(K'+2)0/Tj^ ^^^-[[k’-1)K’0/T^ 


H = (K' + l)e 

(p = 2.88 


( 7 ) 
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for the ground state of nitrogen and 


J = 




( 8 ) 


In reference 5, 9 was introduced as a normalizing factor and is given by 


=I( 

K’ 


^k’+i^p 




(9) 


The number densities in 6 will not have a Maxwellian distribution among the ground 
states but this quantity is still a valid normalizing factor for a particular set of 
conditions. 

At low temperatures there is no appreciable excitation of the vibrational levels of 
other than the V'^ = 0 level so that = N. Therefore, for the 

0-0 band, equation (6) becomes 




FNiq(0,0) 

NH 

9 

%('^R)yU 


which for a given temperature and beam configuration, can be written as 


$(V' = 0,K') = Constant i 



+ J 


where the constant is per second per unit beam current and is independent of 

0 

N since 9 varies as N, By equating the rate of photon emission to the rate of excita- 
tion, Muntz (ref. 5) obtains the emission-time intensity from 


I(K',K’^)^^ = Constant P(K’,K^) 4>(V’ = 0,K’) 


Substituting P^K^Kg) = Pr = Ky(2K’ + 1) into this equation and treating as a con- 
stant gives for the 0-0 band 


' 'exp 


N 


-K 


Qr(T 


K’Ge 




R 


R 


K' 

2K' + 1 


( 10 ) 
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where G is given by equation (2) and Cg is a constant relating the experimental 
intensities to the theoretical intensities for a particular experiment. Define a function 


I(K',K'^) 


= K’Ge 


-K'(K’+1)(#)/Tg 


calc 


Then equation (10) can be written as 


(») 


l{K',iq) 


exp 




N 


K’ 


Q 


r(^ 


R 


I(K’,K") + 

2^calc 2K’ + 1 

/y" 


.1 


(12) 


The reader is referred to Muntz's paper (ref. 5) for more discussion of the steps 
leading to equation (12) since the present derivation is similar. The effect of secondary 
electrons on a particular K' level is introduced through J and depends on the rate of 
secondary excitation to the K" = K' ± 1 levels. The total intensity of the R or 
P branch of the band is not changed by the secondary excitation since one level is 
always populated at the expense of another level. Therefore, the sum over all K' of 
the terms involving J must be zero and from equation (12) 



Let 


and 


Then 




- Cgi 



calc 
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yields 


Normalizing equaticm (12) by 



exp 


”l(K’)1 


"Kk’)] 

_^(0,0) 

exp 

_^(0,0)_ 


calc 




- «lV»)lca.c 


( 13 ) 


where only the upper state quantum number K' denotes the rotational line. 

To evaluate J, values of R^n are required for each ground-state rotational level. 
The selection rule for transitions produced by the low-energy electrons is AK" = ±2 
and R^n is the net value for the transitions between levels K" and K" ± 2. The 
cross sections are denoted by g(v,±2)j^, where the subscript denotes the level from 
which the transition arises and (±2) is AK". 

When the excitation of the K" level from the K" - 2 level is considered, the rate 
of excitation by secondary electrons is 


fR" = '^s^K"-2 r g(v,+2)j^,_2f(v) V dv 


per unit volume per second. To avoid the integral over the unknown f(v), the cross sec- 
tion g(v,+ 2 )j ^,_2 is written in terms of the value for a reference velocity v^ and all 
the cross sections are assumed to vary in a similar manner with velocity. This variation 
does not seem to be unreasonable from the values given in reference 16. Then 


n g(v,+2)^j-n 9 

- "sNk"-2S(''o.+2)K"-2 \ . 

V g(Vo,+2jg.M_2 


Within the accuracy of this assumption, the integral is the same for all K", can be 
treated as a constant, and allows the calculation of relative values for Thus, 

rj^„ = Constant ngN’j^M_2g(voj+2)g^»_2 (14) 


per unit volume per second. 

The rate of excitation required is that due to transitions within the observed gas 
volume. For a field of view of the optics larger than the beam diameter, the beam size 
determines this volume. For large or scattered beams, the field of view of the optics 
usually determines this volume. If p is the beam diameter or field-of-view width, the 
dependence of rj^n on number density is different for the two limits ^ P ^*^*1 
A.g < p. For < p, the excitation rate is given by equation (14) but for Xg > p, the 
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probability of a secondary electron colliding within the observed volume increases as N 
and Tjj" ~ NngN’j^„_2g(vo,+2)jjT._2. The excitation rate Rg,i is the net result of rj^, 

for four transitions and for ^-g > P and is 

Rj^, = Constant Nng |^Njjit^2g(Vo>-2)g;"+2 ^K"-2 ^('^o’‘'’^)k"-2 

- NK”[g(V+2)K" + g('^o^-2)K']) 

Substituting values of R^n into equation (13) and writing the result in terms of K' 
gives for ^g > P 


p(K’)“ 


ri(K’)i 


%,0)_ 

exp 

o 

3 

calc 


^”s^ ^r('^r)vi f{K' + 1)[X(K* + 1)] 


I calc 


(2K’ + 1) 


2K’ + 3 


K'[X(K' - ij 
2K' - 1 


(15) 


and, similarly, for < P 


ri(K’)i 


Fkk’)] 


““s'* 

j 

_^(0,0)_ 

exp 

ho,o)_ 

calc 

”['(0.o3 


V” /(K» + 1)[X(K» + 1)| 


2K’ + 3 


calc 


+ 


K'|X(K’ - l]| 
2K' - 1 


( 16 ) 


where A and B are constants and 

X(K') = N'^._2g(v„,-»2)jj,.2 + N’k.h,2K(V-2)k'*2 

- N'k’[s(''o.+2)K’ + s(v„,-2)k] (17) 

A first approximation to X(K') can be obtained by substituting Maxwellian values 
for N^ti; this substitution gives 
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N 


(18) 


X(K') = 




F(K') 


where 


r _ -(K'-2)(K'-1 )*/tJ 
F(K') = g(v„,+2)^,_^l(2K' -3)e 

[~ -(K’+2)(K’+3)<^/tJ 


(2K' + l)e 


K’+2 

K’(K’+1)0/Tj^ 


p(v+2)g, + g(v-2)g, 


Substituting for X(K') in equations (15) and (16) gives for ^-g > P 


(19) 


'l(K’)' 


ri(K’)i 

_^(0,0) 

exp 

^(0,0)_ 


ilcalc 


An^NK’ 

O 

K’ + 1 


2K’ + 3 



2K’ - 1 


( 20 ) 


and for < p 

O 


ri(K')“ 

a; 

I(K')“ 

_^(0,0)_ 

exp 

^(0,0)_ 


calc 


BngK' 

(K’ + 1)F(K’ + 1) 1 K'F(K’ - 1) 

^(0,0) 

(2K’ + 1)1 

calc 

2K’ + 3 2K’ - 1 _ 


( 21 ) 


For the proposed excitation model, equations (20) and (21) predict that the differ- 
ence between the measured line intensities and values calculated from Muntz's equation 
(ref. 5) varies as NUg at low densities > p) and tends to a dependence on ng at 
high densities (Xg < pj. The transition from one regime to the other occurs for Xg « p 
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and therefore at a density which depends on the geometry of the electron beam and 
optical system. At the present time ng cannot be reliably evaluated as a fxmction of 
the gas density. Estimates can be made of the increase in ng due to the electrostatic 
field set up around the beam (ref. 13), but the quantitative effects of beam spreading and 
the recombination of electrons with ions are unknown. In addition, the constants A 
and B involve the integral over the secondary electron distribution function and cannot 
be reliably evaluated. 

The value of equations (20) and (21) is the predicted dependence of the normalized 
line intensities, and hence the calculated rotational temperature, on the gas density, 
secondary electron number density, and beam geometry. Experimental measurements 
are needed to confirm these predictions and to provide a calibration of the technique for 
a range of temperatures and densities. 

MEASURED ROTATIONAL TEMPERATURES 

The experimental measurements were made in a flowing stream of nitrogen at a 
known temperature using a 30 keV electron beam and beam currents between 0.3 
and 1.5 mA. The impurity level in the nitrogen was 50 parts per million. Figure 2 
shows the heat exchanger used to cool the gas and the position of the electron beam. 



Figure 2.- Schematic diagram of test chamber. 


16 




The heat exchanger is enclosed in a bath of 2 -methyl butane which is thermally insulated 
by the surrounding vacuum. A controlled flow of liquid nitrogen through a coil immersed 
in the 2 -methyl butane permits the cooling of the heat exchai^er to any desired tempera- 
ture between 300 K and 113 K. Alternatively, liquid nitrogen is used in place of a 
2 -methyl butane to obtain temperatures close to 78 K. The heat exchanger consists of 
two chambers separated by a valve which controls the pressure in the first chamber. 

This chamber acts as a reservoir of cooled gas for the main chamber to insure that the 
gas is cooled to the temperatvire of the heat exchanger, even at the lowest pressures. 

The gas temperature is assumed to be the mean value recorded by thermocouples on the 
last baffle and along the exit pipe from the heat exchanger. The diameter of the exit pipe 
was 32 mm and the pumps maintained a maximum velocity through it of approximately 
50 m/sec which is sufficient to prevent convective heating of the gas by the hot drift tube 
and mixing with the warmer gas from the surroimding chamber. (See ref. 8.) 

A calibrated pressure transducer recorded the test chamber pressure through an 
orifice in the chamber wall. The wall temperature remains close to ambient and thermal 
transpiration along the tubing is negligible. The density of the flowing gas is calculated 
from the measured pressure and temperature. 

The electron beam enters the test region through a narrow orifice in the drift tube. 
For pressures below 130 N/m^, the orifice diameter is 0.1 cm but for higher pressures 
is 0.05 cm. Both drift tube and electron gun are insulated from ground potential and a 
grounded shield leaves only a small circular area around the drift tube orifice exposed. 

To measure beam current, the whole test chamber served as the collector and the carbon 
Faraday cup reduced the backscatter of primary electrons from the chamber walls. This 
procedure avoids errors due to beam spreading that arise when only the Faraday cup is 
used as the beam collector. 

Measurements showed that at high beam currents the nitrogen first negative band 
intensities did not vary linearly with beam current, particularly for the higher densities 
and smaller drift tube orifice. The nonlinearity is thought to be due to an interaction 
between the drift tube and the plasma produced by the beam in and around the end of the 
drift tube orifice. The drift tube is slightly positive with respect to ground potential and 
draws a net electron current from the plasma. Unbalanced positive ions left in the plasma 
diffuse to the test chamber walls and thus reduce the indicated electron -beam current 
below the true value. This problem is discussed in reference 8 and leads to errors for 
high beam currents at high densities, but these errors do not seriously affect the conclu- 
sions of this study. 

An optical system formed a 1:1 image of the beam on the entrance slit of a spec- 
trometer which resolved all the lines in the R branch of the 0 -0 band. The slit length, 
oriented perpendicular to the length of the beam image, determined the width of the field 
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of view of the optics. This width was 0.1 cm on either side of the beam center line, 
unless otherwise stated. The slit width determined the resolution along the length of the 
beam so that the volume of gas observed depends on the slit area and the depth of focus 
of the optical system. At low densities, when beam spreading is small, the fluorescence 
is observed from a cylindrical volume defined by the intersection of the field of view Avith 
the beam. For higher densities some of the electrons are scattered beyond the field of 
view (ref. 8). Fluorescence is then observed from a section through the beam of almost 
constant width determined by the slit length. Measurements were made at selected tem- 
peratures for a range of densities by scanning the rotational lines while measuring the 
total fluorescence intensity to monitor for changes in pressure or beam current. 

Figure 3 shows the influence of beam current and density on the measured rota- 
tional temperature at 125 K. The symbols denote the values for a fixed density with 

variable beam current, and Ni correlates the ratio ^ for the limited 

calc/ 

range of beam currents. The solid line is a cubic polynomial fitted through the points by 

the least -squares method. It is noted that (Tp) /T™ does not approach 1.0 as 

"''calc/ 

Ni tends to zero and it increases to an almost constant value for Ni > 6 X 10^® micro- 
amperes X molecules per meter^. 



Figure 3»- Calculated rotational temperature against Ni. ~ 
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The present results for other temperatures are similar to those presented in fig- 
ure 3. The ratio (Tj^) /T^ increases with density and reaches an almost constant 

ca>lc/ 

level at high densities. This level decreases with increasing temperature. 

DISCUSSION 

Equations (20) and (21) predict that AI(K'y^jl^Q pjj ^ varies as ngN at low 
densities > p) and tends to a dependence on ng at high densities ^Xg < p^ . Figure 4 
shows jl(K’)y/l^Q for two rotational levels at 111 K plotted against Ni. At high 

densities |l(K’)yi^Q , and hence AI(K')y^[l(0^o3 ^ , varies linearly with Ni and 

at low densities it varies approximately N times faster, as predicted by equations (20) 
and (21). The transition between the two regimes occurs for a number density of 
approximately 1.5 X 10^2 per meter^, corresponding to a total collision cross section 
of 3 X 10 "19 meter^ which is the same order of magnitude as that given by Takayanagi 
et al. (ref. 16) for electrons of 1 eV energy. Furthermore, 1(3) originating from the 
abundantly populated rotational levels is reduced and 1(13) originating from the 
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sparsely populated levels is increased. Similar trends observed for aU the present data 
strongly support the assumption of low-energy secondary electron excitation of the 
ground-state rotational levels. 

The ratio jl(K'yi^Q does not tend to jl(K’yi^Q as N tends to 

zero for the present measurements. This condition could be due to some unaccounted 
for secondary effect or to the variation of ng with density. Dunn and Self (ref. 17) 
state that for N < 3 X 10 18 molecules per meter^, the beam -generated plasma is 
t3^ically of about the same density as the beam. At higher densities, the electrostatic 
forces become important and references (13) and (15) show that ng is increased by 
several orders of magnitude for N = 3 X 10l9 molecules per meter^. Therefore, Ug 
must increase rapidly with density for 3 x 10l8 < N < 3 x 10l8 molecules per meter^. 

All the present results are for higher densities and would therefore appear to tend to the 
wrong limit when plotted against Ni. 

The potential distribution in the beam -generated plasma governs the distribution 
of ng (refs. 13 and 17) but the magnitude of Ug depends on the rate of production of 
secondary electrons and ions by the primary beam. Although the total beam current does 
not change greatly, the distribution of current density j changes with beam length and 
density because of the spreading of the beam through collisions. Therefore, unless the 
distribution of ng is independent of the distribution of j, the present measurements at 
a distance of 2.8 cm along the beam should reflect the variation of j with density. The 
ratio I(K')y^I^Q varies approximately as N%, tends to Ni at high densities, and 
thus implies that over the observed area of the beam, ng does not depend on the distri- 
bution of j but only on the total rate of production of secondary electrons which is pro- 
portional to Ni. The solution of Bienkowski and Harbour (ref. 13) implies that this 
concept should be approximately true since the secondary electron number density at the 
beam center tends to become independent of beam diameter for large ratios of chamber 
to beam diameter. 

This point is of major importance in any attempt to correct rotational temperature 

measurements by using curves of (Tj,) /T^ against Ni. If ng is not independent 

' ■'calc/ 

of the distribution of j, then the ratio (Tp>) /T^ would vary along the length of the 

' 'calc/ 

beam. The reasonable agreement of the present results (figs. 5 and 6) with those of 
other investigators (refs. 6 and 18) and the evidence of figure 4 suggest that Ug is 
approximately independent of the distribution of j at least for moderate beam lengths. 

Figure 4 shows that there are three density regimes to consider. For 
N < 3 X 10^8 molecules per meter^, a rapid increase in ng with density is thought to 
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Flgiire 5-“ Comparison with Williams' result at j8 K. = 15 • 
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Figure 6.- Comparison with Ashkenas ' results at j8 K using 15 rotational lines. 



be responsible for the large values of AI(K')/Il 


f, , and hence high 
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even at the lowest densities reached in the present experiment. For higher densities, 
ng increases approximately linearly with Ni and the mean free path of the electrons 
varies as 1/N. Therefore, when Xs > p, the number of secondary electrons colliding 
within the observed volume increases as N, and results in AI(K')y|l^Q increasing 

approximately as N%. For Xg < p, all the electrons are observed and AI(K'yi^Q 

increases as Ni. This result explains the variation of (Tjo] jl 

' 'calc/ 


which should be 


a function of Ni when P ^ function of N^i when Xg > p. The correlation 

of the data in figure 3 with Ni for the complete density rai^e is not understood, but the 
range of beam currents is not large. 


When Xg < p, the geometry of the electron beam and optical system are less 
important since practically all the secondary electrons produced within the observed vol- 
ume collide within that volume and (Ttj) /t^ increases only very slowly with Ni. 

'calc/ 

However, when Xg ^ p, the geometry determines how many secondary electrons collide 

within the observed volume and hence the rate at which (Tp^ /T^ approaches the 

' 'calc/ 

almost constant value for Xg < p; therefore, for low densities, increasing the slit length 

should increase (Tt,) /T^. Figure 5 compares (Tj,] /T^ measured at 78 K 

' 'calc/ 'calc/ 

for different spectrometer slit lengths. The limited amount of data and the experimental 

scatter precludes any definite conclusion, but (Td) /T^ appears to increase slightly 

^ 'calc/ 

with slit length. Petrie et al. (ref. 19) observed a similar effect of slit length in a 
Mach 11 flow with static temperatures for 20 K and 35 K. Williams (ref. 18) also made 
measurements at 78 K but in his case the slit dimension, in the direction perpendicular 
to the beam length, appears to be 2 X 10 "5 meter. The experimental scatter is too large 
to draw any conclusion about the effect of the small slit dimension but the results are 
scattered about the present data and increase with Ni. 


Ashkenas (ref. 6) measured rotational temperatures at 78 K and room temperature 

but did not state what beam currents were used. The maximum beam current for his 

data (ref, 6) appears to be 3.5 x 10"^ ampere, whereas for the present results it was 

a maximum of 1.5 X 10 ~ 3 ampere and decreases to about one -fifth of this value at high 

densities. When Xg < p, (Ttj) /t^ is insensitive to Ni and therefore these dif- 

' 'calc/ 

ferences in beam current should not be important. This result is confirmed in figure 6 
which compares the results for 78 K plotted against density only. At low densities 
^Xg > p), the experimental scatter prevents different trends in the two sets of data from 
being detected. 
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In principle, equation (15) or (16) allows the relative values of AI(K')/|I(q 

/I— ^ ^ calc 

to be evaluated, but the secondary-electron excitation redistributes the molecioles among 
the ground-state rotational levels and is not known. A first approximation can be 

obtained by assuming a Maxwellian distribution which leads to equations (20) and (21). 
With equation (19), either equation (20) or (21) gives values of AI(K') /-AI(3) as a func- 
tion of K' only, and these values are compared with experimental values in figures 7 
and 8. Quantitative agreement is not expected because of non -Maxwellian distribution, 
but the qualitative agreement is very good and sparsely populated levels are populated at 
the expense of the more abundantly populated levels. Qualitative agreement at higher 
temperatures is equally good and, as the peak in the rotational distribution moves toward 
higher K", the depopulation of successively higher rotational states is predicted and 
agrees with the data. At temperatures of several hundred K, the low rotational states 
have relatively low populations and secondary electron excitation would increase their 
population at the expense of the higher states. This increase leads to low values of 

(Tp) and could explain Hunter's (ref. 20) results at high temperatures, but the 
'calc 

report does not state which rotational lines were used. 



Fig\ire 7-- Relative intensity error at 79 K- 
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Aerodynamic Measurements 

An important use of the electron-beam fluorescence probe is the study of high 

Mach number low -density aerodsmamics and consideration must be given to the effect of 

high flow velocities. The present results strongly support the hypothesis of excitation 

of ground-state rotational levels by low-energy secondary electrons. The qualitative 

effect of high flow velocities on this model can be deduced by reference to the rarefaction 

parameter 6 defined by Bienkowski and Harbour (ref. 13) as the ratio of the beam 

radius to a modified mean free path. High velocities can only affect the (Tn) /T^ 

' ''calc/ 

through the distribution of Og since the rate of production by primary electrons is inde- 
pendent of velocity. Collisions between fast -moving neutral molecules and ions in the 
beam will sweep the ions downstream and the electrostatic attraction between ions and 
electrons causes a redistribution of the secondary electrons. Bienkowski and Harbour 
(ref. 13) concluded that for 6 = 1, there should be little effect of velocity, but for 

6 » 1, the distribution of ng, and hence (Tt, ) affected. A value of 6 of 1 

' ''calc/ 

corresponds to approximately 10 N/m2 at room temperature for many electron-beam 
situations and most of the present results are in the regime where 6^1. 

Comparison of the present results with wind-tunnel measurements for known flow 
conditions and similar beam parameters would determine the effect of flow velocity. 
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Measurements are available in a Mach 6 flow (ref. 21) for a static temperature of 71.5 K 

and 6 of order 1. Extrapolating the present results to this temperature gives 

(Ttj) = 1.13 ± 0.03 compared with 1.09 ± 0.03 for the wind-tunnel data; thus, 

' 'calc/ 


it is confirmed that there is little effect of the velocity for 6 « 1. For all other known 
wind-tunnel measurements, there are either geometric differences or the beam currents 
are not given. In addition, there is some doubt about Ni as a correlating parameter. 

If a beam current of 4.5 X 10 ampere is assumed, the Mach 9 data of Williams, 
Hornkohl, and Lewis (ref. 22) agree with the present data to within 6 percent. Their 
Mach 10 data are considerably higher, possibly because of uncertainties in the flow con- 
ditions. The Mach 3.92 data of Robben and Talbot (ref. 9) for a beam current of 
10 “2 ampere are very much lower than the present results for the same Ni but show 
reasonable agreement when compared only in terms of the density. Maguire's free-jet 
measurements at 154 K (ref. 7) cover a range of densities for 6 > 1, but the beam cur- 
rents were not given and the slit was parallel to the beam image. The effect of the dif- 
ferent geometric arrangement is less important at high densities, and for 
N = 1.23 X 1023 molecules per meterS^ Maguire's result is 4 percent lower than that of 
the present study. This difference shows that there is no strong effect of flow velocity 
on 


(Tj,) . Such a result is to be expected since (Tt,) /' 

^ ^^calc ^^calc/ 


T^ is insensitive to ng 


at high densities (kg « p] and therefore to the influence of flow velocity on the distribu- 
tion of ng. For lower densities, the discrepancies increase to 10 percent for 
N = 1.32 X 1 q 22 molecules per meter^, but no conclusion can be drawn because of the 
different geometries. 


Correction Procedure for 

Many variables influence (T^) and simple correction procedures will apply 

^'calc 

only for specific conditions. The present measurements used beam currents of approxi- 
mately 10 "2 ampere and a field of view of the optical system of 0.1 cm on either side of 
the beam center line. The drift tube orifice was 0.1 cm or 0.05-cm diameter, but the 
beam spreads beyond this at the measurement station which is 2.8 cm from the drift tube. 

A correction procedure using curves of (Tj^^ /T^ as a function of Ni will apply 

calc/ 

only for beam currents for which the correlating parameter Ni is valid. The results 
might be different for beam lengths considerably greater than 2.8 cm, and at low den- 
sities (kg > p^, the field of view is important. 

For each set of data at a fixed temperature T^, a least -squares cubic polynomial 

was fitted through the values of (Td) /T^ as a function of Ni. The polynomial 

' 'calc/ 


(Tr) 


calc 
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TABLE I.- POLYNOMIAL COEFFICIENTS FOR 


Coefficients for T of 
I 167 K 210 K I 



79 K 

111 K 

125 K 

1 

167 K 

210 K 




K 

■max “ 


Co 

1.226 

1.103 

1.110 

1.047 

1.028 

Cj X lo21 

15.11 

9.86 

I 2.24 

j 

2.66 

4.20 

C 2 X 1053 

‘ -22.69 

-16.0 

1.30 

0.174 

-4.71 

C 3 X 10'^9 

11.65 

8.85 

-2.37 

1 -2.07 

1.95 


Cg X 1079 1 12.74 


^max 


1.091 

1.106 

1.037 

1,025 



8.84 

0.597 

3.40 

4.46 

2.26 

7.7 

-15.3 

4.82 

-2.15 

-5.36 



9.23 

-4.64 

-0,224 

2.30 


18,58 








Cq 

1.162 

1.105 




1.063 


Cj X 102^ 

12.81 

6.14 

1,45 

3.51 

4,64 

1.77 

5.23 

Cg X 1033 

-24.06 

1 

-9.46 

2.35 

-1.44 

-5.73 

-0.373 

-11.20 

Cg X 10'79 

14.45 

5.41 

-2.89 

-1.19 

1 

2.51 

-0.623 

i 

8.78 


Cj X 102'7 
Cg X 10^^ 
Cg X lo”^® 


Cj X 1Q2'7 
Cg X 1053 
Cg X 10'79 


^max 


1.142 

1,126 

1.103 

1,020 

0.992 

1.054 

11.81 

4.00 

0.345 

3.78 

7.82 

1,66 

-23.07 

-6.03 

0.438 

-2,46 

-13,23 

0.898 

14.06 

3,88 

-4.09 

-0.217 

7.58 

-2.27 


^max 


(Ni)max ^ ^ 10^'’ 

(Ni)^.„ 7.8 X 1023 


Range 

7.7 X 1025 
1.5 X 10^4 


1.033 

5.13 

-14.05 

13.55 


1.121 

1.126 

1.075 

0.995 

0.978 

1.042 

0,996 

11.77 

2.83 

2.54 

5.11 

7.67 

2.22 

8.81 

-24.34 

-2.07 

-0.581 

-4.82 

-12.7 

2.65 

-24,97 

15.37 

0.401 

-0.954 

1.26 

7.55 

-5.69 

23.40 


j^. microamperes X molecules 
meter^ 

9 X 10^5 7.2 X 10^5 8 X 10^5 7.2 X 10^3 
2 X 10^4 2 X 10^^ 2 X 10^4 1 X 10^^ 


1.02 

10.22 

-34.28 

34.77 

1.00 

11.04 

-34.06 

32.22 


6 X 1025 6 X 102E 

1.3 X 1024 6.4 X 10 


coefficients are given in table I, where is the number of spectral lines used 

except for 300 K, when the first two lines were neglected because of overlapping by 
the P branch. Within the range of Ni given at the bottom of this table for each tem- 
perature, values of (Tj,^ /T^ can be calculated from 

'■ 'calc/ 

(Tr)^^1c/Tw = % + ClNl C2(N«2 „ C3(N1)3 

Values of (Tj^) /t^ calculated from these polynomials are plotted in figure 9 for 

calc/ 

selected values of Ni. 



The calculated rotational temperature depends on the density, but to measure den- 
sity when quenching is appreciable, the temperature must be known and an iterative pro- 
cedure is required to correct the measurements. (See ref. 8.) First estimates are made 

of and Nj and equation (1) used to calculate and 

calc calc/ 

computed from the polynomials by substituting Nj^i. There is some advantage in norma- 


lizing (Tjj) /t^ by the value at 300 K, that is, 
'calc/ 



) 

300 K 


since the 


normalized ratios are insensitive to small changes in Ni, A polynomial is fitted through 
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the normalized temperature ratios as a function of temperature and stored on the com- 
puter for use in successive iterative steps. If R(T) denotes the normalized temperatur 

ratio at temperature T^, then [(Tp) 1 '^'^ ~^('^i)[(^r) 

L calc/ J,p^ |_ calc/ J 300 K 

obtained by substituting Tj and Nji into the relevant polynomials. A better estimate 
of the temperature is Tj^ = fixed, the loop is 

repeated by substituting this improved estimate into the polynomials until the calculated 
temperature converges to a constant value T 2 . A second estimate of the density N 2 
can then be obtained by using T 2 and the quenching cross sections given in reference 8. 
Repeating this process gives successive estimates of temperature and density. If the 
estimate Nj was reasonable, the original polynomial for R(T) can be used for each 
iterative step without making serious errors. This procedure involves more computation 
than that given previously (ref. 8) but allows a greater range of Ni to be covered at low 
temperatures. 

Other correction procedures were investigated and one consequence of the proposed 
secondary-electron excitation model is that ^ AI(K'^jl^Q = 0. This summation 

varies slowly with the error in the temperature. Combined with large errors in AI(K’), 
this condition makes it impractical to correct the temperatures by minimizing the sum- 
mation. A simple and apparently successful procedure developed earlier (ref. 10) was 
based on the incorrect assumption of direct excitation of the states by secondary 

electrons. Subsequent study shows that the equations in reference 10 are approximate 
functions fitted through Ashkenas’ data. (See ref. 6.) The successful application of this 
procedure in reference 21 is the result of similar beam currents and geometry for the 
two experiments. Correcting this shock-layer data (ref, 21) by the present method did 
not change the temperatures sufficiently to affect the conclusions of that paper. 

CONCLUDING REMARKS 


Rotational temperatures measured by the electron-beam fluorescence technique 
increase with density and beam current and depend on which spectral lines are used. 

The secondary mechanism responsible for the high rotational temperatures cannot be 
conclusively identified, but excitation of the ground-state rotational levels by low-energy 
secondary electrons can account for all the observed results. This excitation mechanism 
redistributes the molecules among the ground-state rotational levels and thus increases 
the population of the sparsely populated levels at the expense of the more abundantly 
populated levels. 


28 


For this excitation model, the spatial distribution and number density of the sec- 
ondary electrons and the field of view of the optical system determine the rotational- 
temperature error for a given density. High flow velocities can influence the measure- 
ments through collisions with ions, and thus result in a different spatial distribution of 
electrons and ions around the beam. For low densities, when the secondary-electron 
mean free path is large compared with the beam diameter or field-of-view width, there 
is little influence of high flow velocities. At high densities, when the secondary -electron 
mean free path is relatively small, the calculated temperature is insensitive to the 
secondary-electron number density, and hence to all but strong effects of the flow on the 
distribution of the secondary electrons in the observed gas volume. For these two 
regimes the present correction scheme should be reliable, but at intermediate densities 
the flow will reduce the secondary electron concentration in the observed gas volume 
and the temperature will be underestimated. 

Further studies are needed to establish reliable correlating parameters for the 
calculated temperatures and to determine more accurately the effect of high flow veloc- 
ities, particularly at intermediate densities. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 26, 1971. 
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